Signature of entangled eigenstates in the magnetic response of two coupled flux qubits. 
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We study dissipative dynamics and a magnetic response of two coupled flux qubits interacting 



with a high quality tank in the framework of the impedance measurement technique (IMT). It is 
shown that the observation of the difference between a sum of IMT signals from separated qubits 
and the signal from the system when both qubits are in the degeneracy point (IMT deficit) implies 
D ■ immediately a formation of the entangled two-qubit eigenstates. 

Q 

The impedance measurement technique (IMT) 1 ' 2 has proven to be among effective methods for experimental inves- 
CN . tigation of flux quantum bits. In the framework of this approach the qubit is inductively coupled to the high quality 
tank (LC circuit) that modifies the impedance of the circuit and, therefore, the angle between the tank voltage and 
the bias current applied to the tank. This modification is proportional to the second derivative of the qubit's energy 
or to the magnetic suceptibility of the qubit 3 . The IMT method was succesfully employed for studying the single 
i-Sh . flux qubits in equilibrium 4,5 and strongly non-equilibrium states 6 . Recently, the IMT measurements were performed 
for the system of two inductively coupled flux qubits 7 with independent control of the bias in each qubit. The areas 
of the qubits were practically the same. Because of this, it was necessary to apply an additional flux created by a 
supplementary wire to change a relative position of the IMT dips on a graph describing a current-voltage dependence 
as a function of a common bias produced by the dc component of the current in the tank. Applying the additional 
flux allows to bring the IMT dips from the different qubits into coincidence (overlapped regime) as well as to pull 
them apart (regime of separated qubits). Effects of interaction and entanglement of the qubits have been observed by 
comparing IMT signals from the separated qubits with the response of the two-qubit system when both qubits were 
in the degeneracy point. The IMT dip for the system in the overlapped regime is shown to be smaller than the sum 
of IMT amplitudes from two separated qubits (IMT deficit) that points to the formation of entangled eigenstates of 
two coupled qubits. 

Here we present detailed calculations of the magnetic response of the two-qubit system coupled to a dissipative 
environment. Besides that we elucidate a relation between the IMT deficit and an entanglement of pure eigenstates 
of the system. 

We consider two coupled flux qubits interacting with a dissipative environment and with a high-quality LC circuit 
(tank) . The Hamiltonian of a such system has the form H = Hq + Ht + Hint + Hdiss , where the Hamiltonian of the 
\Q , two-qubit system is 

CN 

H = -A o0 £°> - A b a^ + e a a^ + e h of + J^of , (1) 

co ; 

Ht is the Hamiltonian of the tank with the current It, the qubit-tank interaction is given by the term 

23 '. H ™t = -(Aa4 0) + V-W)J T , ( 2 ) 

and H diss = -ai a 'Q a ~ af 
and Q b . The IMT signal 3 ' 7 , 



and Hdiss = —&z Qa — &zQb describes the weak coupling of the qubits to the dissipative baths with variables Q a 



tane = -2^ X >T), (3) 



is determined by the real part of the magnetic susceptibility of the double-qubit system, x'(u;), where 

In the general case the eigenstate \fi),fi = 1,2,3,4, of the Hamiltonian H : H \fi) = E u \/j,), can be represented as a 
superposition of the double-qubit states in the standard basis 

= ci„|00> + c 2 J01) + c 3M |10) + c 4m |11), (5) 

where |00) = |0) a ® |0) b is a direct product of states of the qubits "a" and "b" with cri a) |0) a = - |0} a , of 1 |l) a = |l) a , 
and so on. The states \fi) and \v) are orthonormal: = 5^ v . 
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We introduce matrix elements for and : 

*V = (»\(ri a) W) = -c*n c iv - c *2^ c 2u + c* 3 ^c 3 „ + c* 4lM c 4v , 

*V = (t-i\(ri b) \v) = -<?xyC\ v + c^c 2 „ - c^c 3u + c* 4ll c iu . (6) 

With these notations the matrices and af^ can be represented as 4 x 4 matrices: af^ 1 — u^ v p^ v \ — 

Y^iiv v iivP^ ,iV \ where p(^) = \p){v\ is the matrix with matrix elements p^jp — 5 ail 5p v . This matrix has only one 
non-zero (unit) element located on the /w-place. Then, the Hamiltonian H can be represented as 

H = J2 E aaP {aa) - ]T(Qa + XaI T )u a0 p^ + J^iQb + A 6 /t) W Q/3) - (7) 

a af3 a/3 

The matrices p^ v ^ can be considered as Heisenberg operators which obey the equations: 

if>W = -CJ^p^ - £>„ Q Q„ + V va Q b ) P M + £>a M Qa + V ali Q b ) P W - 
a a 

+ \ b V va )I T p^ + ]T(Aa?W + >*Vc,»)lTP (aV) , (8) 
a a 

with uj^ u = E^—E v . For the response of the heat baths on the action of the qubits we obtain from the theory of open 
quantum systems 3,8 

Qa(t) = Q { °\t) + J dt l {u a , tj ,ip aa {tM) + v a >f3>iPab{tM)}p {a ' P ' ) (h)}, 

Qb(t) = Q ( b 0) (t) + J dhiu^pipbafah) + v a , p ,ip bb (t,t 1 )}p^'P'\t 1 )}. (9) 
Here , are unperturbed variables of the heat baths, 



Vaa(t,tl) = <i[Qi 0) (t).Qi 0) (*l)]->e(t-tl)= / ^ 



e-^-^Xoaiu) (10) 
is the linear response function and the susceptibility of the "a" heat bath, 

*i) = (i[Q(°>(t),Qi 0) (ti)]_>(9(t-ti) (11) 

is a response function describing a cross-correlation between "a" and "b" heat baths with the corresponding sus- 
ceptibility x a h(w), 6{t) is the Heaviside step function. We introduce also corresponding correlation functions and 
spectra 

M aa (t,h) = ((l/2)[Q^(t),Q^( tl )} + )= J ge-^'-'^M, 
M ab (t,h) = <(l/2)[Q<°>(t),Q$ 0) (ti)] + > = J ^e-W-VS^w) (12) 
together with combinations 

<Pu'f(t,h) = (iKvQ^it) + v va Qf\t), v^i 0) (ti) + v a , ,Qf\u)]-)6{t - h), 

M v'a' (Mi) = <(l/2)[u„ Q Q(°>(t) + v„ a Q< b °\t), u a ^Q<°>(ti) + v a ,p,Qf\t 1 )] + ). (13) 

With these definitions the "collision" terms in Eq. (8), averaged over the thermodynamically equilibrium state of the 
heat baths, take the form 

((l/2)[u va Q a (t)+v va Q b (t),p^ a \t)]+) = 
J dt 1 M^\t,t 1 )(i[p^(t),p^^(t 1 )}.) + J di 1 </(t,i 1 )((l/2)[p(^)(i),p (a ' /5 ' ) (ii)] + ); 

((l/2)[u aft Q a (t)+v aft Q b (t),p av (t)} + ) = 
J dt.MZfitMmP^it),^'^^)]-) + J dt 1 ^'(t,ti)<(l/2)^(t),p( a '' 3 ')(ti)] + ), (14) 
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where M(r) = M(t)6(t). In the case of weak interactions between the qubits and the baths we can use a free evolution 
of the matrices p^ a \t) = e^*"^*"* 1 ) p^ a ) (t x ) and p^(t) = e ^-(*-*i )p(«") (£ x ) to calculate (anti)commutators in 
the collision terms (14) using the relations: 

[p (^) (tl)) pto'u') (ti)]± = (5q , q/9 (^') (tl) ± ^,^(«'«) ( tl ) 



[p^(tl),p (a ' /3 ' ) (il)]±=5 a VP W ' , (il)±VaP^^(il), 

where the all Heisenberg matrices are taken at the same moment of time t\. As a result we obtain the following 
equation for the averaged operators (p^^ ) (we consider thereafter the evolution of the averaged variables only, and, 
because of this, omit the brackets (..}): 



(«/3')/ 



4. a ' v )i 



lU^p 



+ i 



a 

a/3 



It 



E 

a/3 

E 

a/3 

E 

a/3 



dti ( M - -if 



{t-t 1 )e w ^ {t - tl) p w) {t 1 ) + 



Pp. 



dh ( M + -ip ) (t- tiie^V-^pW^ih) + 



dh [M - -Lp 



VOL 

up 

CtfJ, 



{t - h)e w ^ {t ^ p {afi \h) - 



E 

ap 



dh[M+ -<p 



Pa 



(t - ^) e ^-(t-ti) p (M( tl ) 



(15) 



with (M — ^<p)"f = — ^Vva- in the equilibrium state, when 7t = 0, the averaged matrices p(^) = p^ 1 ^ with 



/i^i/ are equal to zero, whereas for the diagonal matrices /?q we obtain the equation 



E 



s + -x) Ka)+ [s--x\ {0J av ) 



(ota>) 

Po = 



/T } E 



(w„ Q ) + [S + -x) (u av ) 



(16) 

For the heat baths with the same parity regarding the time inversion we have: Xab{u) — Xba{<jj), S a b(^) = Sba{w), 
where the spectrum S a b{u>) is related to the imaginary part of the corresponding susceptibility Xabi^) according to 
the fluctuation-dissipation theorem: S a b = x'abi 10 ) coth(w/2T). Here T is the temperature of the heat baths "a" and 
"b". The same relations take place for S aa (ui) and x"„(w), as well as for Sbb(w) and x'bb( w )- I n view of the fact that 
u av = u* a ,v av — v* a and with the notation p 1 ^'^ = p v the equation (16) can be rewritten as 



E{^ 



coth 



~2T~ 



- 1 



coth 



{\u va \ 2 X'La(u va ) + KalM&^a) + (u va V* a + VuaKJXabi^ua)} = 0. 



(17) 



The solution of this equation is given by the Gibbs distribution: p v = cxp(—E u /T)/Z, where Z = J2a ex P(~E a /T). 
The term, related to the interaction of qubits with the tank's current It in Eq.(15), transforms to the form: 
i-i^aUvfj, + ^bV V n){Pn — Pv)It- The interaction with the tank drops out from the equation for diagonal elements 
p^ vl/ \ so that {5p( uv \t)/5lT{t')) = 0. It means that only non-diagonal elements, p^ v \ give non-zero contributions to 
the susceptibility x( w ) (4)- For the Fourier transform of the functional derivative (8p^ v \t) / 'dixit')) we obtain from 
Eq. (15): 



0» 



'lu + Ef, - E v + iT^iuj) 



ipu - Pix)- 



(18) 



A small frequency shift of the double-qubit system caused by the dissipative environment is omitted here. With a 
notation 



k va {u) = \u va \ Xaai") + W a \ Xbbi^) + iUvcxVta + V„a<Jx"bM 



(19) 
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the frequency-dependent decoherence rate is given by the expression 



E 1 

^ 2 



H,.l,|^| - I 



A„ a (u; +w M c) + 

A„ a (uj + u au ). (20) 



Because of the frequency dependence of T^ v {lj) the dissipative evolution of the averaged elements, (p^ v \t)), is 
governed by the decoherence rate, T^ U (E V — E^), which differs from the decoherence rate V ^ v {lot) incorporatED 

into the low-frequency susceptibility x(i^t) (4) of the coupled qubits. Taking into account the relations — 

u llv p( llv } — J2^ v v^p^ v \ for the magnetic susceptibility we find the equation 

{Ku^ + A b ?v)( — — (w)} = 

OlT 

{XaU^ + XbV^iXaU^ + XbV^). (21) 



E 



, UJ + En - E v + lYnv (to) 



If we recall that and v^v are the matrix elements of the and (6) between the eigenstates and \v) of 
the Hamiltonian Ho (1) and neglect the resonant frequency of the tank compared to the energy differences between 
the qubit's levels, lot *C E v — E^, then the current-voltage angle in the tank (3) is given by the expression 

tan6 = -4^ (22) 

where 

+A a A b ( M |a( a ) \v){y\af \p) + X a X b (p\a^\y) (i/|a<»> (23) 

For the weak qubit-bath couplings the decoherence rates, r ^(ujt) <C E v — E^, drop out of Eq.(22), which describes 
the IMT signal. 

The first line in Eq. (23) corresponds to the contribution of individual qubits to the IMT signal. It is clear that this 
contribution survives even when the eigenstates \v) are products of single-qubits states. The terms in the second 
line of Eq. (23) are related to the simultaneous flipping of both qubits. This term vanishes when all eigenstates of the 
two-qubit system are non-entangled. To show that we suppose that the states \p), \v) can be represented as products 
of the wave fuctions related to the qubits "a" and "b": \p) = \p) a \p)b, \v) = \v) a \ v )b- It is the case, for example, for 
the separated qubits, when one qubit is in the degeneracy point, but another - far out of degeneracy. The states \p) 
and \v) are orthonormal, so that for p ^ v we have (p\v) = (p\v) a (p\v)b = 0. It means that {fj\v) a = or (p\f)b = 0, 
or both. Then the matrix elements {^u^ \v), {v\<j^\p) involved in the expression for the IMT deficit 7 , 

(iVW = \ a \ b {(M a) W)(v\4 b) \v) + (n\ai b) W)(v\4 a) \»)}, (24) 

have the form: (p,\a^\v) — (p,\a^\v) a (p,\v)b and (i/\ai b ^\p) = (v\p) a^a^ \p)b- The products of these matrix elements 
contain the term {n\v)b{v\y) a which should be equal to zero due to the normalization condition, {p\v) a {ii\v)b = 0. 
We can conclude, therefore, that non-entangled eignestates of the two-qubit system gives no contribution to the IMT 
deficit {Rfj, v )def- Thus, the observation of non-zero IMT deficit directly points to a formation of entangled eigenstates 
in the system of two coupled flux qubits. It should be noted, however, that zero IMT deficit does not mean that all 
eigenstates of the system are non-entangled. 

In summary we have performed calculations of the magnetic susceptibility and IMT signals for two coupled flux 
qubits. We have shown that the observation of non-zero IMT deficit presents a sufficient condition for the existence 
of entangled eigenstates in the two-qubit system. 

Discussions with E. Il'ichev, A. Izmalkov, M. Grajcar, M. Amin, A. Maassen van den Brink, and A. Zagoskin are 
gratefully acknowledged. I especially thank Alexandre Zagoskin for critical reading of the manuscript. 
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